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Abstract 

We study the properties of the field in the fundamental mode HEn of a vacuum- 
clad subwavelength-diameter optical fiber using the exact solutions of Maxwell's 
equations. We obtain simple analytical expressions for the total intensity of the 
electric field. We discuss the origin of the deviations of the exact fundamental mode 
HEn from the approximate mode LPoi- We show that the thin thickness of the fiber 
and the high contrast between the refractive indices of the silica core and the vacuum 
clad substantially modify the intensity distributions and the polarization properties 
of the field and its components, especially in the vicinity of the fiber surface. One 
of the promising applications of the field around the subwavelength-diameter fiber 
is trapping and guiding of atoms by the optical force of the evanescent field. 
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1 Introduction 



Optical fibers with diameters from micrometers to millimeters have found 
many important applications in technology. The guiding properties of such 

* Corresponding author. Tel.: +81-424-43-5476; Fax: +81-424-43-5507. 

Email address: fam@kiji.pc.uec.ac.jp (Fam Le Kien). 
1 Also at Institute of Physics and Electronics, Vietnamese Academy of Science and 
Technology, Hanoi, Vietnam. 



Preprint submitted to Elsevier Science 



1 February 2008 



fibers have been well studied [1]. Most fibers use core and cladding materials 
with nearly the same refractive indices. Such fibers are referred to as weakly 
guiding fibers. The description of weakly guiding fibers is greatly simplified 
due to the use of the linearly polarized (LP) modes, which are approximate 
solutions of Maxwell's equations and are superpositions of two corresponding 
nearly degenerate eigenmodes. 

Recently, thin waveguides have become attractive for a wide range of potential 
practical applications. It has been pointed out that the intense evanescent 
field in the vicinity of a tapered fiber can be used as an atomic mirror [2]. 
Generation of light with a supercontinuum spectrum in a thin tapered fiber 
has been demonstrated [3]. The evanescent waves from zero- mode metal-clad 
subwavelength-diameter waveguides have been used for optical observations of 
single- molecule dynamics [4]. Several types of dielectric submicro meter- and 
nanometer-diameter wire waveguides have been fabricated, and their guiding 
properties have been investigated [5]. It has been proposed to use the optical 
force of an evanescent wave around a thin fiber for atom trapping and guiding 
[6]. Thin fiber structures can be used as building blocks in the future micro- 
and nano-photonic devices. 

Unlike the case of optical fibers with diameter larger than the wavelength, 
many properties of subwavelength-diameter fibers have not been adequately 
investigated. In such a thin fiber, the original silica core is vanishing. Therefore, 
the original silica clad acts like a core while the surrounding vacuum (or air) 
acts like a clad. Due to the high refractive-index contrast between the silica 
core and the vacuum clad, the description of subwavelength-diameter fibers in 
terms of weakly guided LP modes is questionable. It has been demonstrated 
that such fibers have interesting properties such as high power density at the 
fiber surface and cylindrical asymmetry in the field distribution [2], large pen- 
etration length of the evanescent wave [6] , enhancement of the power fraction 
of the field outside the fiber [6,7], and large waveguide dispersion [7]. How- 
ever, the intensity distributions and polarization orientations of the fields in 
subwavelength-diameter fibers have not been studied. 

In this paper we study the field intensity distributions and polarization orienta- 
tions in the fundamental mode HE n of a vacuum-clad subwavelength-diameter 
optical fiber using the exact solutions of Maxwell's equations. We show that 
the thin thickness of the fiber and the high contrast between the refractive in- 
dices of the silica core and the vacuum clad substantially modify the intensity 
distributions and the polarization orientations of the field and its components. 
The study of such properties is necessary not only for academic interest but 
also for practical applications such as trapping of atoms by the optical force 
of an evanescent wave around a thin fiber [6] . 

Before we proceed, we note that, in related problems, the decay of an atom 
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in the presence of a fiber was considered for the first time by Katsenelenbaum 
in f949 [8]. Excitation of fiber modes with a dipole source and influence of 
a fiber on the decay rate of a single atom were investigated in [8,9,10,11]. 
In the present paper, we do not consider how the modes can be excited. 
Due to the adiabatic tapping condition [12], the fundamental mode HEn 
of a subwavelength-diameter fiber can be excited by the standard coupling 
techniques [5]. 

The paper is organized as follows. In Sec. 2 we describe the general model of a 
thin fiber and examine a fundamental mode with quasi-linear polarization. In 
Sec. 3 we study a fundamental mode with rotating (circulating) polarization. 
Our conclusions are given in Sec. 4. 



2 Thin fiber and fundamental mode with quasi-linear polarization 



Consider a thin single-mode optical fiber that has a cylindrical silica core of 
radius a and refractive index n\ and an infinite vacuum clad of refractive in- 
dex n 2 = 1. Such a fiber can be prepared using taper fiber technology. The 
essence of the technology is to heat and pull a single-mode optical fiber to 
a very thin thickness maintaining the taper condition to keep adiabatically 
the single-mode condition [3,12]. Due to tapering, the original core is almost 
vanishing. Therefore, the refractive indices that determine the guiding prop- 
erties of the tapered fiber are the refractive index of the original silica clad 
and the refractive index of the surrounding vacuum. The refractive index and 
the radius of the tapered silica clad will be henceforth referred to simply as 
the fiber refractive index rt\ and the fiber radius a, respectively. 

We send a light of wavelength A, frequency u, and free-space wave number 
k = 2n/\ = uj/c through the fiber. Under the condition V = ka\Jn 2 — w% < 
V c = 2.405, the fiber can support only one mode, referred to as the funda- 
mental mode HE n . The longitudinal propagation constant (3 of this mode is 
determined by the eigenvalue equation [1] 
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Here the parameters h = [n\k 2 — f3 2 ) 1 ^ 2 and q = {(3 2 — n 2 k 2 Y^ 2 characterize 
the fields inside and outside the fiber. The notation J n and K n stand for the 
Bessel functions of the first kind and the modified Bessel functions of the 
second kind, respectively. 
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We first study a fundamental mode with quasi-linear polarization. In the cylin- 
drical coordinates (r, tp, z), the solutions of Maxwell's equations for the Carte- 
sian components of the electric field E in such a mode are given, for r < a 
(inside the fiber), by [1] 



E x = -^[(1 - s)J (hr) cosv?o - (1 + s)J 2 (hr) cos(2y9 - <p )] e l ^ z \ 

E y = - s)J (hr) sin v^o - (1 + s)J 2 (hr) sin(2y9 - <p )] e { ^ z \ 

E z = A.h(hr) cos(v? - <p ) e*^"^ , (2) 

and, for r > a (outside the fiber), by [1] 



,_, . 3 JAha) r/ . T ^ , . 

E x = -iA— (1 - s)K (qr) cos (p 

2qK 1 (qa) 



+ (1 + s)fT 2 (gr) cos(2y9 - <p )] e^'"^, 

E ^4 q W) [(1 - s)Ko{qr)s ^° 

+ (1 + s)K 2 (qr) sin(2y9 - </? )] e^"^, 
E z = A^j^-K^qr) cos(y? - <p ) e^'\ (3) 

where s = [(qa)~ 2 + {ha)~ 2 ]/[J' l {ha) / haJi(ha) + K[(qa) / 'qaKi(qa)]. The coeffi- 
cient A is determined by the normalization condition. The angle (po determines 
the orientation axis of the polarization of the field. The two sets of solutions 
corresponding to <po = and <p Q = n/2 express two different polarizations, 
aligned along the x and y axes, respectively. 

For practical applications such as trapping of atoms by the optical force of an 
evanescent wave around a thin fiber [6], it is necessary to calculate the optical 
potential, which is proportional to the total intensity \E\ 2 of the electric field. 
A rigorous expression for \E\ 2 can be easily obtained with the help of Eqs. (2) 
and (3). For r < a, we find 



\E\ 2 = 9inUo(hr) + uJ 2 (hr) + fj 2 (hr) 

+ [uJ 2 (hr) - f p J (hr)J 2 (hr)] cos[2(y9 - <p )]}. (4) 

For r > a, we get 

\E\ 2 = g out {K 2 (qr) + wK 2 (qr) + fK 2 2 (qr) 

+ [wKl{qr) + f p K (qr)K 2 (qr)} cos[2(y9 - <p )]}. (5) 



4 



Fig. 1. Total intensity \E\ 2 and x-component intensity \E X \ 2 of the electric field 
in an exact fundamental mode HEn with a quasi-linear polarization, along the 
x direction, for the parameters of a conventional weakly guiding fiber. The two 
profiles are indistinguishable because the intensities of the y and z components of 
the field are negligible. The parameters used: a = 4 /im, A = 1.3 /im, n\ = 1.4469, 
n 2 = 1.4419, and (p = 0. 

Here we have introduced the parameters u = 2h 2 /j3 2 (l — s) 2 , w = 2q 2 /(3 2 (l — 
s) 2 , f = (l + s) 2 /(l — s) 2 , and fp = 2(l+s)/(l — s). We have also introduced the 
notation g in = \A\ 2 /2u and g ont = \A\ 2 J 2 (ha)/2wKf(qa). The terms Jl(hr) 
and Kl(qr) in the expressions (4) and (5), respectively, correspond to the total 
intensity of the electric field in the mode LPoi. The other terms describe the 
deviations of the exact mode HEn from the approximate mode LPoi. 

The above expressions are mathematically valid for any core radius a and 
any pair of refractive indices n\ and n 2 where n\ > n 2 . Before applying these 
expressions to thin tapered fibers, we recall the case of conventional single- 
mode fibers, where A = (ni — n 2 )/ni "C 1 and a > A. In this case, we have [1] 
s = -1 + (qa/V) 2 [haJ (ha) / 'Ji(ka)]A + 0(A 2 ) = -1. Consequently, the in- 
dependent terms containing + J 2 {hr) and (1 + s)K 2 (qr) in the expressions 
(2) and (3) are negligible. Hence, the transverse component E y or E x will be 
zero if ipo = or ip = n/2, respectively. In addition, we have h,q -C /3, so the 
longitudinal component E z is small. Thus, an exact HEn mode with quasi- 
linear polarization of a conventional weakly guiding fiber can be approximated 
by a LPoi mode [1]. The electric field in the LPoi mode is given by E = 
Fee^-M, where F = AJ (hr) for r < a and F = A[J (ha) / K (qa)}K (qr) 
for r > a. The polarization vector e can be chosen as e = x or e = y. The 
intensity distribution of the LP i mode is cylindrically symmetric. 

In Fig. 1, we plot the total intensity \E\ 2 and the x-component intensity \E X \ 2 
of the electric field in an exact fundamental mode HEn with a quasi-linear po- 
larization, along the x direction, for the parameters of a conventional weakly 
guiding fiber. As seen, the transverse profiles of \E\ 2 and \E X \ 2 are indistin- 
guishable. Thus the polarization of the field is almost completely linear. We 
also observe that the two profiles are almost perfectly cylindrically symmetric. 

We now study the case of vacuum-clad subwavelength- diameter fibers, where 
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the relations A 1 and a > A are not satisfied. In this case, the factor f + s 
is not negligible. More importantly, the decay parameter qa of the evanes- 
cent wave may become sufficiently small that K^qa) and K 2 (qa) are much 
larger than K (qa). Hence, the terms containing K^qr) and K 2 (qr) in the 
expressions (3) may become significant in the outer vicinity of the fiber sur- 
face. Due to the factor K\(qr), the longitudinal component E z outside the 
fiber may become substantial. Due to the factor K 2 (qr), the terms containing 
the trigonometric functions cos(2(/? — ip ) and sin(2(^ — ip ) in the expressions 
(3) may become significant. These terms lead to the azimuthal dependences 
of the transverse components (E x , E y ) of the field outside the fiber. In ad- 
dition, when qa is small and n\ is much different from n 2 , the parameter h 
may become comparable to (or even larger than) the parameter (3. Due to 
this fact, the component E z of the field inside the fiber may also be not neg- 
ligible compared to the components (E x ,E y ), respectively, see Eqs. (2). Thus 
the longitudinal component E z may be substantial in both regions r < a and 
r > a. So we see that the properties of the exact fundamental mode HEn may 
become substantially different from the linearly polarized mode LP i. 

To demonstrate the features of vacuum-clad subwavelength- diameter fibers, 
we perform numerical calculations. The single-mode condition for vacuum- 
clad silica-core fibers requires that a/A < 0.36, that is, the fiber radius be 
about three times smaller than the light wavelength. We consider the case 
where the fiber radius a is small compared not only to the light wavelength 
but also to the evanescent-wave penetration length A — l/q. Such thin fibers 
can be used to trap atoms. Indeed, it has been shown that the optical force 
of a red-detuned evanescent wave around an optical fiber can balance the 
centrifugal force under the condition qa = a/A < 0.93 and that this condition 
can be achieved for vacuum-clad silica-core fibers with a/A < 0.28 [6]. We 
study such thin fibers in the rest of this paper. 

To be specific, we choose the fiber radius a = 0.2 jum and the light wavelength 
A = 1.3 fim for simulations. The corresponding refractive indices of the silica 
core and the vacuum clad are ri\ = 1.4469 and n 2 — 1, respectively. The solu- 
tion of the exact eigenvalue equation (1) yields ha = 1.0075, qa = 0.0827, and 
Pa = 0.9702. The normalized size parameter is V = ka^Jn\ — n\ = 1.011 < 
2.405, indicating that the considered fiber is a single-mode fiber. For the above 
parameters, we find s = —0.9937. Since the fiber radius a is thin compared 
to the light wavelength A, the penetration length of the evanescent wave is 
A = 12 a, large compared to the fiber radius [6]. Because of this, a majority 
of the field intensity distribution is in the outside of the fiber [6,7]. 

In Fig. 2, we plot the total intensity \E\ 2 of the electric field in a fundamental 
mode with quasi-linear polarization. We choose the x axis as the major ori- 
entation axis of polarization (ip = 0). As seen, the behaviors of the profiles 
inside and outside the fiber are very different from each other. A conspicuous 
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Fig. 2. Cross-section profile of the total intensity \E\ 2 of the electric field in a 
fundamental mode with quasi-linear polarization. The inset shows the inner part 
of the profile, which corresponds to the field inside the fiber. The parameters used: 
a = 0.2 /im, A = 1.3 //m, m = 1.4469, ni = 1, and <po = 0. 

discontinuity of the field intensity is observed at the fiber surface. This discon- 
tinuity is due to the boundary condition for the normal (radial) component 
E r of the electric field. The high contrast between the refractive indices n\ 
and ri2 of the silica core and the vacuum clad, respectively, makes this effect 
dramatic. 

As seen from Fig. 2, the spatial distribution of the field intensity is not cylin- 
drically symmetric at all. A strong dependence of the field intensity on the az- 
imuthal angle is observed in the outer vicinity of the fiber surface. This depen- 
dence is due to the terms wKf(qr) cos[2({p— ip )] and f p K (qr)K 2 (qr) cos[2(ip — 
ipo)] in Eq. (5). These terms are substantial because the small magnitude of 
the parameter qa = 0.0827 makes K\(qr) and ^(gr) dominant over Ko(qr) 
in the vicinity of the fiber surface. Indeed, we have Ki(qa) / Ko(qa) = 4.6 
and K2(qa) / Ko(qa) = 111.7. In the inner vicinity of the fiber surface, the 
azimuthal dependence of the field intensity is moderate. 

To get a look at the spatial distributions of the field along different radial 
directions, we replot in Fig. 3 the intensity \E\ 2 of Fig. 2 as a function of x 
at y = (as a function of r at ip = 0) and as a function of y at x = (as a 
function of r at ip = vr/2). For comparison, we plot by the dashed lines the 
intensity |-Elp| 2 of the corresponding approximate mode LP i, which is given 
by |-Elp| 2 = 9mJo(hr) for r < a and |-E LP | 2 = g out K^(qr) for r > a. The 
comparison between the region x/a > 1 of Fig. 3(a) and the region y/a > 1 of 
Fig. 3(b) shows that, outside the fiber, \E\ 2 decays along the x direction faster 
than along the y direction. We observe that the radial dependence of \E\ 2 along 
the y direction has a small peak in the outer vicinity of the fiber surface. Far 
away from the fiber surface, the decays of \E\ 2 along the x and y directions tend 
to have the same behavior as that of the intensity |-Elp| 2 of the approximate 
mode LPoi- The comparison between the region x/a < 1 of Fig. 3(a) and the 
region y/a < 1 of Fig. 3(b) shows that, inside the fiber, \E\ 2 decreases along 
the x direction slower than along the y direction. A discontinuity of \E\ 2 is 
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Fig. 3. Total intensity \E\ 2 of the electric field in a fundamental mode with 
quasi-linear polarization as a function of x at y = (a) and as a function of y 
at x = (b). For comparison, the intensity of the electric field in the corresponding 
approximate mode LPoi is shown by the dashed lines. The parameters for this figure 
are the same as those for Fig. 2. 

observed at the fiber surface in the x direction, but not in the y direction. 

To get insight into the field components of the fundamental mode, we plot 
in Fig. 4 the cross-section profiles of the intensities |-E X | 2 , \E y \ 2 , and \E Z \ 2 
of the Cartesian-coordinate components of the electric field. In addition, we 
plot in Fig. 5 these intensities as functions of the azimuthal angle (p. The 
comparison between the scales of the vertical axes of Figs. 4(a), 4(b), and 
4(c) shows that the intensity \E y \ 2 of the minor transverse component E y 
and the intensity \E Z \ 2 of the longitudinal component E z are weaker than 
the intensity \E X \ 2 of the major transverse component E x . However, \E y \ 2 and 
| E z | 2 are not negligible at all. Therefore, the alignment of the total electric field 
vector E = E x x + E y y + E z z may be substantially deviated from the major 
orientation axis x of polarization. According to Eqs. (2) and (3), the major and 
minor transverse components E x and E y of the electric field have the same 
phase. Therefore, the vector orientation of the total transverse component 
Ej_ = E x ~k + E y y does not vary in time. Thus the transverse component Ej_ of 
the field is linearly polarized with respect to the time evolution at each fixed 
spatial point (x,y). Meanwhile, the phase of the longitudinal component E z 
differs from that of the transverse components E x and E y by tt/2. Due to this 
phase difference, the total electric field E rotates elliptically with time, in a 
plane parallel to the fiber axis z. 

The inset in Fig. 4(a) and the dotted curve in Fig. 5(a) show that, inside the 
fiber, \E X \ 2 practically does not vary with (p. The inset in Fig. 4(b) and the 
dotted curve in Fig. 5(b) show that, inside the fiber, \E y \ 2 is negligibly small, 
so the orientation of the vector E^ almost coincides with the x axis. Thus the 
total transverse component E^ is not only linearly polarized with respect to 
the time evolution at each fixed local point but also is almost linearly polarized 
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Fig. 4. Cross-section profiles of the intensities |.Ke| 2 j l-E^I 2 ) an d \E Z \ 2 of the 
Cartesian-coordinate components of the electric field in a fundamental mode with 
quasi- linear polarization. The insets show the inner parts of the profiles, which cor- 
respond to the field inside the fiber. The parameters for this figure are the same as 
those for Fig. 2. 

in the space inside the fiber. 

The main profiles in Figs. 4(a) and 4(b) as well as the dashed and solid curves 
in Figs. 5(a) and 5(b) show that, outside the fiber, \E X \ 2 and \E y \ 2 substantially 
depend on ip. This behavior is quite different from the behavior of the field 
inside the fiber. 



Figures 4(c) and 5(c) show that the longitudinal-component intensity \E Z \ 2 
is substantial not only in the outer vicinity of the fiber surface but also in 
the inner vicinity, unlike the minor-transverse-component intensity \E y \ 2 . In 
addition, \E Z \ 2 substantially varies with ip not only outside but also inside the 
fiber, unlike the major-transverse-component intensity l-E^ 2 . The azimuthal 
dependence of \E Z \ 2 is profound in both sides of the fiber surface because E z 
is proportional to cos(y2 — <po), see Eqs. (2) and (3). 
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Fig. 5. Azimuthal profiles of the intensities l-E^I 2 , l-Eyl 2 , and \E Z \ 2 of the Carte- 
sian-coordinate components of the electric field in a fundamental mode with 
quasi-linear polarization. The parameters for this figure are the same as those for 
Fig. 2. 
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Fig. 6. Orientation angle 8 = arctan[Re (E y )/~Rje {E x )\ of the transverse component 
Ej_ of the electric field in a fundamental mode with quasi-linear polarization. The 
parameters for this figure are the same as those for Fig. 2. 

At each fixed spatial point, the vector orientation of the total transverse com- 
ponent Ej_ of the electric field does not, as discussed earlier, vary in time. 
However, the orientation of Ej_ may vary in space, especially in the outer 
vicinity of the fiber surface, where the minor transverse component E y is sub- 
stantial. We plot in Fig. 6 the orientation angle 9 = arctan[Re (E y )/Re (E x )] of 
Ej_ as a function of (p. The dotted curve in Fig. 6 shows that, inside the fiber, 
the orientation angle 9 almost does not vary with <p. This feature confirms that 
the total transverse component Ej_ inside the fiber is almost linearly polarized 
in space. It is due to the fact that the (^-dependent terms in the expressions 
for the transverse components inside the fiber are negligible. Meanwhile, the 
dashed and solid curves in Fig. 6 show that, outside the fiber, the orientation 
angle 9 does vary in space. For r/a = 1.5, the maximal azimuthal variation of 
9 is about 0.06 n (see the dashed curve). The comparison between the solid 
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and dashed curves in Fig. 6 shows that, when r increases from a, the azimuthal 
variation of 9 reduces. These features are due to the fact that the ratio between 
the factors K 2 (qr) and K (qr) of the ^-dependent and -independent terms is 
large for qr <C 1 but monotonically reduces to unity with increasing r. 



3 Fundamental mode with rotating polarization 



We now study a fundamental mode with rotating (circulating) polarization. 
In the cylindrical coordinates, the solutions of Maxwell's equations for the 
cylindrical components of the electric field E in such a mode are given by [1] 



E r = F r e ±ilp e i{u>t - pz \ 
E {p = ±F {p e ±i * j^-M, 

E z = F z e ±iv j^-M. (6) 

Here the functions Fj (j = r, ip, z) describe the radial dependences of the field 
components. They are defined, for r < a, as 

F r = -^[(1 - s)J (hr) - (1 + s)J 2 (hr)}, 
F v = A^[(l - s)J (hr) + (1 + s)J 2 (hr)}, 

F z = AJ 1 (hr), (7) 
and, for r > a, as 



Fr = ~ lA 2qY^) [{l ~ S)K ° iqr) + + S ^ qr ^ 
F * = A YqK^) [{1 ~ S)K ° iqr) + S ^r)], 

The upper (lower) sign in Eqs. (6) corresponds to the clockwise (counterclock- 
wise) circulation of photons around the z axis. According to Eqs. (6), we have 
I Fj | = \Fj\. According to Eqs. (7) and (8), the functions Fj are independent of 
the azimuthal angle tp. Hence, the intensities \Ej\ 2 of the cylindrical-coordinate 
components of the field are independent of ip, and so is the total intensity \E\ 2 
of the electric field. 
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According to Eqs. (7) and (8), we have F r /\F r \ = —iF^/lF^. Then, it follows 
from Eqs. (6) that E r /\E r \ = ^fiE v /\E v \, that is, the two complex quadra- 
tures E r and E v have a difference of n/2 between their phases. Therefore, the 
polarization of the total transverse component Ej_ = E r r + E^ip of the field 
is either circular or elliptical. With increasing time, the real part of the com- 
plex amplitude vector Ej_ rotates along a circle or an ellipse. The semimajor 
and semiminor axes of the ellipse are aligned along the unit vectors f and (p, 
respectively, and are equal to \E r \ = \F r \ and \E V \ = \F V \, respectively. When 
\E r \ = \E V \, the elliptical polarization of the transverse component of the field 
becomes almost circular. 

We can easily calculate the total intensity \E\ 2 of the electric field in a fun- 
damental mode with rotating polarization. For the field inside the fiber, we 
obtain 

\E\ 2 = 2g- m [J 2 (hr) + u.J 2 (hr) + f J 2 2 (hr)}. (9) 

For the field outside the fiber, we get 

\E\ 2 = 2g out [K 2 (qr) + wK 2 (qr) + fK 2 2 (qr)}. (10) 

The terms Jl{hr) and Kl(qr) in the expressions (9) and (10), respectively, 
correspond to the total intensity of the electric field in the mode LPoi. The 
other terms describe the deviations of the exact fundamental mode HEn with 
rotating polarization from the approximate mode LP i. The comparison be- 
tween Eqs. (4) and (9) and between Eqs. (5) and (10) shows that the total 
intensity of the electric field in a fundamental mode with rotating polariza- 
tion is the sum of the corresponding intensities for two constituent modes 
with quasi-linear polarizations. The (/^-dependent terms cancel each other and 
therefore do not appear in Eqs. (9) and (10). 

The above expressions are mathematically valid for the fundamental mode 
with rotating polarization of a fiber with an arbitrary core radius a and an ar- 
bitrary pair of refractive indices rii > n 2 . To demonstrate the features of 
vacuum-clad subwavelength- diameter fibers, we perform numerical calcula- 
tions for the parameters of the previous section, namely, for a = 0.2 /im, 
A = 1.3 /an, rii = 1.4469, and n 2 = 1. 

In Fig. 7, we plot the cross-section profile of the total intensity \E\ 2 of the 
electric field in a fundamental mode with rotating polarization. As seen, |.E7| 2 is 
azimuthally independent, that is, is cylindrically symmetric. The discontinuity 
of the field at the fiber surface, created by the boundary condition and the 
high contrast between ri\ and n 2 , divides the profile into two parts. The inner 
part of the profile is hidden behind the outer part and is shown separately in 
the inset of the figure. 
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Fig. 7. Cross-section profile of the total intensity \E\ 2 of the electric field in a 
fundamental mode with rotating polarization. The inset shows the inner part of 
the profile, which corresponds to the field inside the fiber. The parameters for this 
figure are the same as those for Fig. 2. 
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Fig. 8. Radial dependence of the total intensity |-E| 2 of the electric field in a funda- 
mental mode with rotating polarization (solid line). For comparison, the intensity 
|£lp| 2 i obtained for the corresponding approximate mode LPoi, is plotted by the 
dashed line. The parameters for this figure are the same as those for Fig. 2. 

In Fig. 8, we plot by the solid line the total intensity \E\ 2 of the electric field 
as a function of the radial distance r. For comparison, we plot by the dashed 
line the intensity | .Exp 1 2 of the corresponding approximate mode LP i, which 
is given by |-Elp| 2 = 2fl , in^o(^ r ) for r < a and |-Elp| 2 = 2g ou t-Ko (l r ) 101 r > a. 
As seen, the difference between \E\ 2 and |-Elp| 2 is small but not negligible in 
the region 0.3 < r/a < 3. Inside the fiber, \E\ 2 (solid curve) decreases slower 
than |-Elp| 2 (dashed curve). Such a behavior of |-E"| 2 is due to the contributions 
of the last two terms in Eq. (9). These terms contain the functions J\ and J2, 
which increase in the region of small argument. In the outer vicinity of the 
fiber surface, \E\ 2 (solid curve) decays faster than |-Elp| 2 (dashed curve). Such 
a behavior of \E\ 2 is due to the contributions of the additional terms wK 2 (qr) 
and fK 2 (qr) in Eq. (10), which decay faster than the basic term K^(qr) in the 
region of small argument. When r/a is large enough (r/a > 3), the difference 
between \E\ 2 and |-Elp| 2 can be neglected. 

The intensities |-Er-| 2 , l-E^) 2 , and |.Ez| 2 of the cylindrical-coordinate components 
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Fig. 9. Intensities l-EVl 2 , |-£v| 2 , and |-Ez| 2 of the cylindrical-coordinate components 
of the field in a fundamental mode with rotating polarization. The parameters for 
this figure are the same as those for Fig. 2. 

of the field are determined by the cylindrically symmetric functions \F r \ 2 , \F V \ 2 , 
and \F Z \ 2 , respectively. We plot these intensities as functions of r in Fig. 9. As 
seen, the field intensity distributions inside (r/a < 1) and outside (r/a > 1) 
the fiber have very different behaviors. Due to the boundary condition and the 
high contrast between the refractive indices of the silica core and the vacuum 
clad, the normal (radial) component E r has a conspicuous discontinuity at the 
fiber surface. The tangential components E v and E z are, however, continuous. 

Unlike the intensities of the cylindrical-coordinate components E r , E v , and 
E z , the intensities of the Cartesian-coordinate transverse components E x and 
E y are, in general, not cylindrically symmetric. In Fig. 10, we plot the cross- 
section profiles of the intensities l-E^I 2 , \E y \ 2 , and \E Z \ 2 . In addition, we plot 
these intensities in Fig. 11 as functions of (p. 

The insets in Figs. 10(a) and 10(b) as well as the dotted curves in Figs. 11(a) 
and 11(b) show that, inside the fiber, \E X \ 2 and \E y \ 2 are almost equal to each 
other and practically do not vary with ip. The reason is that the 99-dependent 
terms in the expressions for F r and F v in Eqs. (7) are negligible. Because 
of this, we have F r = —iF v and, consequently, E r = ^iE^. Hence, we find 
E x S ir re <M-^) ) E y 2* ±iF r e i ^ t ~ f3z \ and consequently E x S ^%E y . Thus the 
two orthogonal quadratures E r and E v as well as E x and E y have almost equal 
magnitudes and a relative phase difference of tt/2. This indicates that, inside 
the fiber, the total transverse component Ej_ of the field is almost circularly 
polarized. 

The main profiles in Figs. 10(a) and 10(b) as well as the dashed and solid 
curves in Figs. 11(a) and 11(b) show that, outside the fiber, the intensities 
of the transverse components E x and E y are different from each other and 
substantially vary with <p. The comparison between the dashed lines and the 
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solid lines in Figs. 11(a) and 11(b) shows that the azimuthal dependences of 
\E X \ 2 and \E y \ 2 outside the fiber reduce with increasing r. 

The terms ±(1 + s)K 2 (qr) in Eqs. (8) lead to \F r \ ^ \F^\ and, hence, to 
\E r \ ^ \Ep\. Due to the substantial difference between \E r \ and \E V \ in the 
outer vicinity of the fiber surface, the polarization of the transverse component 
of the field in this region is truly elliptical. We recall that the axes of 
the elliptical circulation orbit are aligned along the unit vectors r and (p. 
Consequently, the orientation of the orbit varies with increasing ip. We also 
recall that the lengths of the axes of the ellipse are equal to \E r \ = \F r \ and 
l-E^I = \F V \. Hence, the ellipticity of the orbit is proportional to the factor 
(1 + s)K2(qr)/(l — s)Ko(qr). This factor reduces with increasing r. Thus both 
the orientation and the ellipticity of the orbit of polarization circulation vary 
in space. This behavior is different from the case of conventional light beams 
with elliptical or circular polarization. Furthermore, we note that the ellipticity 
of the orbit of polarization circulation does not depend on ip. This indicates 
that the orbit rotates circularly in space. 

Figures 10(c) and 11(c) show that the longitudinal-component intensity \E Z \ 2 
is perfectly cylindrically symmetric, as expected, in the whole cross-section 
plane. Although \E Z \ 2 is small compared to \E X \ 2 and \E y \ 2 , it is not negligible 
in the vicinity of the fiber surface. 



4 Conclusions 

We have studied the properties of the field in the fundamental mode HE n of 
a vacuum-clad subwavelength- diameter optical fiber using the exact solutions 
of Maxwell's equations. We have obtained simple analytical expressions for 
the total intensity of the electric field and have identified the origin of the 
deviations of the exact fundamental mode HEn from the approximate mode 
LP i. We have shown that the thin thickness of the fiber and the high contrast 
between the refractive indices of the silica core and the vacuum clad substan- 
tially modify the intensity distributions and the polarization properties of the 
field and its components, especially in the vicinity of the fiber surface. 

We have examined the case of a fundamental mode with quasi-linear polar- 
ization. We have shown that a substantial azimuthal dependence of the total 
intensity is observed in the vicinity of the fiber surface. The transverse compo- 
nent of the field is linearly polarized in time at each fixed local point. However, 
the total electric field vector rotates elliptically with time, in a plane parallel 
to the fiber axis. Inside the fiber, the transverse component is not only linearly 
polarized in time but also almost linearly polarized in space. Outside the fiber, 
the orientation angle of the transverse component of the field varies in space. 
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We have also studied the case of a fundamental mode with rotating (circu- 
lating) polarization. We have shown that the total intensity is azimuthally 
independent, that is, is cylindrically symmetric. We have found that the total 
intensity of the electric field outside the fiber decays faster than that of the 
approximate mode LP i. The difference between the exact and approximate 
modes is relatively small in the case of rotating polarization. Consequently, 
the underlying physics of the optical potential of the evanescent wave around 
a vacuum-clad subwavelength-diameter fiber is basically the same as that of 
the approximate mode LP 01 . However, the magnitude of the difference is not 
negligible in the vicinity of the fiber surface. Therefore, the use of the exact 
solutions of Maxwell's equations is required in a systematic quantitative treat- 
ment for a thin fiber. We have shown that, inside the fiber, the total transverse 
component of the field is almost circularly polarized. In the outer vicinity of 
the fiber surface, the polarization of the transverse component of the field is 
elliptical. In this region, the orientation of the orbit of polarization circulation 
rotates in space while the ellipticity of the orbit reduces with increasing ra- 
dial distance. This is different from the case of conventional light beams with 
elliptical or circular polarization. Our results are helpful for studying and 
developing new miniaturized high-performance photonic devices. One of the 
promising applications of the field around the subwavelength-diameter fiber is 
trapping and guiding of neutral atoms by the optical force of the evanescent 
field. 
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